TXTLE OF THE INVENTION Qf 

Method and system for detecting and classifying 
analog and digital tdecorronunicauons sig^s- 

The invention relates m geu telecom munications signals, 

typically: CW, FM. PSK, AM, DSB-SC, BPS*, 
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BACKGROUW. OF THE WVWIKW 
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Other, more recent, onor art « ^as, each yielding a binary 

clas5ifiMti on decisions - • ^ such methods * » a decision 

output. As described to refer— (21. W. - d "* 
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performance accuracy. intolerance to carrier frequency 
publish** to-date utherent d ^ ^ 

errors, their erratic performance at low » and men 
) schemes reliably under real-time operations. : 

Other forms of classification techniques are described m the folio mg J? 
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OTHERS by using the whole vector of six binary demons as an rnput to a logt ^ 
OTHERS, by g computation of the vector 

The drawback to this method is that it usually pen 
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— ^^trrpUability distribution of the input signal ampli.de * 
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^ r^R cmals The combined information is 
„ AM one to FM steals -d one o ^ 

— to . pre—d se^s — ^Z* ^ - P— 

, v^c-o ccd psK. The main difference oeiweeu 
AM, FM, DSB, SSB, K>*- m ~ egsi ka of the feature vector, as 

invention is the computationally-intense parallel process 

* t ™*1 nrocessin- of the vector which is less computer hungry- 
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centre frequency ottneoverausp 

, spared to those of typical modulation schemed such as AM. FSK, PSK or SSB. 
) compared _ representation of the signal 

main drawback of this method is that it uses only P 

to perform its computation. 

In referee U 0]. a method is nsed to discriminate bn- an I'M signal and a 
W4.DQPSK in analo g AMPS and distal DAMPS sy.cn.- The variation u, 
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signals that are encountered in High Definition 
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fa reference [13]. the specW energy dtstribution 0 f<he input sign* U compared 
wiU, the pre-stored energy distributions of FDM/FM signals containing speetfic 
parameters. Recognition of . specific form of si^l is declared if the inpu, dtsmbutron 
notches one of the pre-stored versions. T* method is limited to . sing.e form of srgnal 
fearore and to a very specific modulation format. 

In reference [14], Ore demodulated signal of an FM receiver is Cashed accordmg 
,o the voice coding algorithm »* Pressed it nd. method differs from the preset 
invention because edification is applied on the demodulated stgnal. 

Jn Terence [151, a method is used to generate a decision-tree cl^fier from a set 
of records. It differs from the present invention in that It does no, consider me specrfic 

classification of modulation formats. 

Reference [16] describe a method and apparatus for detecting and classing 
signals ttiat arc the additive combination of a few constant-a^pHtude sinusal 
components. The main drawback of this scheme is that it cannot be appUed to the 
modulations treated under the present invention, except for CW. 

In reference [17], a sequence of estimated magnitudes is generated fium the 

received signal at me symbol ^ . 

„„u™A Coital data modem signals. This method is limited to a 
representation of known voiceb and digital data moaem s 

, , ^e»rit«rhp series of binary decisions performed under tlte 

single decision level, as opposed to trie series 01 omdi, 
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f~ U SUMMARY OF THE INVENTION % 

f It is now an object of the invention to provide a more straightforward and 

1:3 computationally simpler method of modulation recognition than neural network based 
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20 methods. 

IS It is another object of the invention to provide a modulation recognition method 

1 5 which directly exploits the fundamental principles of the decision tree based methods in 

° order to automatically and accurately perform recognition of a wide variety of modulation 

formats that are embedded in unknown communications signals. These modulation 
formats comprise the following set: CW. AM, FM, FSK, DSB-SC, BPSK, QPSK, u/4- 

QPSK, MPSK, NOISE, and OTHERS. 

It is a further object of the present invention to provide a modulation recognition 
method that successfully classifies the signal's embedded modulation amid SNRs as low 
as 5 dB, carrier frequency errors up to +/- 50% of the sampling rate, and carrier phase 
errors up to +/- 180 degrees - all without a priori knowledge of the symbol transition 
timing of the signal. 
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Tt is a further object of the invention to extract the digital complex baseband of 
unknown signals that are measured off-air in real-time, on-line in real-time, or from 
storage, and then accurately classify the signal's embedded modulation through an orderly 

series of decision functions- 

It is still a further object of the invention to provide a system for recognizing the 
type of modulation of a modulated signal having a carrier frequency comprising: 

(a) a receiver section for extracting from the modulated signal a complex baseband 
signal; 

pi (b) a pre-classification stage for generating a pre-processed signal from the 

baseband signal; 

(c) means for examining amplitude variations in the pre-processed signal to 

identify the type of envelope thereof as being one of a constant envelope and 
one of an irregular envelope; 
0 1 (d) means for estimating the carrier frequency error in the pre-processed signal; 
L (e) means for correcting the pre-processed signal for the carrier frequency errors to 

generate a carrier-corrected signal; 
P / (f) means for classifying the type of modulation in the carrier-corrected signal, 

based on the type of envelope identified in the examining step. 
0 \ (g) means for classifying the type of modulation embedded in a constant envelope 
signal according to the phase and frequency contents of the pre-processed and 
carrier-corrected signals, such modulations being categorized within the set 
{CW, FSKand FM}; and 
Q I (h) means for classifying the type of modulation embedded in an irregular 

envelope signal according to the phase and amplitude contents of the carrier- 
corrected signal, such modulations being categorized within the set {AM, 
DSB-SC, BPSK, QPSK, n/4-QPSK, MPSK, OTHER}. 

f The method of the invention provides a unique decision tree architecture that 
► automatically performs recognition of a wide variety of embedded modulation formats in 
unknown communications signals. The principal method of the invention extracts the 
digital complex baseband of the unknown signal and then determines and classifies, 
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through an orderly series of signal processing factions, the signal's embedded 
modulation scheme with a high degree of accuracy. Without knowledge of the symbol 
transition timing, the method of the invention can be used for successfully performing 
modulation recognition of measured signals amid SNR's as low as 5 dB, carrier frequency 

5 errors of up to +A50% of the sampling rate, and carrier phase errors of up to +/- 1 80 
degrees. The tolerance to carrier frequency errors is limited only by the frequency 
bandwidth of the unknown signal, since the only requirement is to prevent the frequency 
error from shifting the input signal outside of the sampling bandwidth of ± half of the 
sampling rate. The preferred method according to the invention is computationally simple, 

1 0 making it possible to observe the signal during a time frame of preferably less than 1 00 
msecs. The modular architecture of the method of the invention allows for expansion of 
new modulation classifiers. 

1 BRIEF DESCRIPTION OF THE DRAWINGS 

1 5 Exemplary embodiments of the invention will now be described in more detail 

with reference to the appended figures, in which the referenced numerals designate similar 

/ 

parts throughout the figures thereof, and wherein: 

Figra-e J.'is a top-level block diagram that illustrates an embodiment of the present 
invention for niea^fnng and classifying an off-air unknown signal. 

Figure illustrates, in a top-level flow chart, the modulation classification section 

* shown in Figure y 

Figure ^illustrates a sample of the spectral power density of the signal and 
embedded noise that can be measured in the pre-cJassification stage shown in Figure 1 arid 
detailed in Figured. 

25 Figure Estates in a graph the degradation in the classification success rate of 

four specific signal types versus the normalized carrier frequency errors, without any 
provision for correcting the frequency errors in accordance with mis invention. 

FigureSui a block diagram illustrating the signal processing functions in the error 
correction block 26 of Figure 2 for correcting the residual carrier frequency error inherent 

30 in a constant envelope signal. 

Figure /is a flowchart illustrating the signal processing functions in the classifier 
block 30 of Figure 2 for classifying the modulation of constant envelope signals. 



20 



8 



t 



Figure ^isa block diagram illustrating the signal processing functions for 
generating the instantaneous frequency components of a constant envelope signal in the 
phase processor block 42 of Figure 6. 

Figure S^sa block diagram illustrating the signal processing functions in the error 
5 correction block 27 of Figure 2 for correcting the residual earner frequency error mherent 
in some irregular envelope signals. 

Figure^isa graph illustrating the performance of the signal processing used in 
Figures 5 and 8 for correcting the residual carrier frequency error inherent in some 
selected signals. ^ 

10 Figure \p is a flowchart illustrating the signal processing functions used for 

classifying the modulation of irregular envelope signals in the classifier block 31 of 
Figure 2. 

Figure U^is a summary flowchart depicting the entire modulation classification 

i 

!J process of the present invention, which entails the combination of figures 2, 6 and 10. 

, 15 / ! 

Qt CL ujc DETAILED DESCRIPTION OF THE INVENTION 

i;n /is 

$ Glossary of Acronyms ! 

I j f For con venience, a glossary of acronyms used in the description of the present 
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invention is given below: 

• 




AM 


Amplitude Modulation 




AMPS 


Advanced Mobile Phone Services 




ASK 


Amplitude Shift Keying 
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AWGN 


Additive "White Gaussian Noise 




BPSK 


Binary Phase Shift Keying 




BT 


Product of a filter bandwidth B and the symbol period T 




BW 


Bandwidth 




CW 


Continuous Wave 




DSB 


Double Sideband 




DSB-SC 


Double Sideband Suppressed Carrier 




FFT 


Fast Fourier Transform 
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FSK 

FM 

GMSK 

GSM 

IS-54 



MPSK 

MUX 

PSD 

PSK 

QAM 

QPSK 

ti/4-QPSK 

RMS 

SNR 

SSB 

VCO 

VSB 



Fluency Shift Keyin< 
Frequency Modulation 
Gaussian Minimum Shift Keying 

Global System for Mobile (Pan European Digital Cellular) 
Interim Standard #54 of die Telecommunications Industry 
Association (also known as US Digital Cellular) 
M-ary PSK 

Multi -Channel or Multiplexing 
Power Spectral Density 
Phase Shift Keying 
Quadrature Amplitude Modulation 

Quaternary Phase Shift Keying 

n/4-shifted Quaternary Phase Shift Keying 

Root Mean Square 

Signal to Noise Ratio 

Single Side Band 

Voltage Control Oscillator 

Vestigial Sideband 



Terminology 

Throughout the description of the present invention the following terms are used in 
accordance with their respective definitions given below: 

Linear vector quantization "Vector quantization" is a process by which a vector, 

defined over a continuous multi-dimensional vector space, is 
assigned one of a limited number of pre-determined 
representative vectors. Quantization occurs at the vector 
level because an infinite number of vector possibilities is 
mapped onto an ensemble containing a finite number of 
vectors. The mapping is performed by dividing the original 
vector space into a number of regions, each of which 
corresponding to one of the finite number of pre-determined 
representative vectors. "Learning vector quantization" is the 
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process that performs the adaptive segmentation of the 
infinite vector space in the predetermined regions daring a 
training period or during normal operation of Hie quantizer. 



Training set The training set is an ensemble of known vectors, signals, ox 

other inputs provided to a system that estimates certain 
parameters over the inputs. The system's internal state is 
thus allowed to adapt to the training set so that it can 
thereafter perform parametric estimations on unknown 
inputs. 

Centroid frequency The centroid frequency of a given power spectral density is 

the frequency corresponding to the "centre of mass" of such 
power spectra) density. It is a purely geometric property of a 
given density and is expressed as: 

^ j / 1 0 where the c,s are the Discrete Fourier Transform (DFT) 

coefficients of the input signal and the /s are the 
corresponding frequencies. 



Zero padding In the computation of the Discrete Fourier Transform (DFT) 

of a given discrete-time signal vector x(n) of finite length N, 
zero padding refers to the appending (the padding) of a 
certain number L of samples (all equal to zero) at the end of 
the vector x(n)- The DFT is then computed on a new 
discrete-time signal vector [x(n) 0 0 0 ...0], of length N-rL. 
This process has the effect of increasing the frequency 
resolution of the original N-point DFT, since it is then 
computed for N+L different frequencies over the same 
frequency range. 
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Unwrapped phase Unwrapping the phase of a complex signal vector refers to 

the process of making the corresponding modulo -2ti phase 
vector continuous over the range from minus infinity to plus 
infinity, by adding multiples of + or -2tt as needed. 

Digital complex baseband The complex baseband representation of signals (also 

referred to as "complex envelope representation") serves the 
purpose of mathematically expressing a real bandpass signal 
in a format which contains the amplitude and phase 
information of the signal but not its carrier or reference 
frequency. This representation is an extension of the familiar 
two-dimensional phasor (vector) representation of sinusoidal 
signals. 
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Consider a real bandpass signal with a narrow bandwidth 
concentrated around a reference or carrier frequency / c - This 
real signal can be mathematically expressed as: 




s(0 = a(t)cos[27f c t + <HO] 

where aft) is the signal amplitude and (pft) is the time^ 
varying phase. The modulation of the signal is contained in 
its amplitude and its phase (the frequency modulation is 
equal to the time derivative of the phase). The signals aft) 
and $ft) show all the modulation informadon carried by &ft), 
and must be preserved in the complex baseband 
representation. The above equation can also be expressed as: 

s(t) = a(t) cosfcKOlcos P^ e r] - a(t) *ntf(0] 
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Substituting: 
L~ v(0 = a(r)sin[ «HOJ- 



we then have a signal: 



/ / O which is the complex baseband equivalent of s ft) and 

contains all the amplitude and phase information of the real 
bandpass signal. 



13 



Overview of Present Invention 

Figure 1 is a top-level block diagram that shows one aspect of the present 
invention - the measurement of off-air unknown signals. This figure includes an antenna 
1, a receiver section 2 and a modulation classification section 3. The receiver section 
includes an amplifier and filter 4, a frequency down-converter 5, a signal sampling 
processor 6 5 and an analog lo digital (A/D) converter 7. Each of these receiver circuits are 
generally recognized in the prior art. The unknown signal, which may comprise a plurality 
of channels, is collected off-air at the antenna 1 and then amplified and filtered by step 4. 
This function also serves to isolate the channel containing the desired carrier frequency 
(f c )> The filtered signal is then down-converted in frequency in step 5 to the analog version 
of the complex baseband signal. While this process reduces f Q to near zero Hz, residual 
carrier frequency errors still remain within the converted signal and must be dealt with 
prior to the classification process in order to ensure that the modulation format is 
accurately recognized. From the down-converter 5, the anaJog complex baseband signal is 
sampled by step 6 at regular time intervals of T s = l/F 5 seconds, and the resulting time 
samples are quantized by the A/D converter 7 to produce a series of quantized time 
samples known as the digital complex baseband signal. Quantization is performed within a 
finite number of levels corresponding to the number of bits of the A/D converter. The 
digital complex baseband signal 8 is then forwarded to the modulation classification 
section 3 ? the principal embodiment of this invention, where the digital signal processing 
for classifying the modulation schemes begins. 

The modulation classifier 3 of Figure I illustrates a simplified view of the 
modulation classification process, consisting of a pre-classification stage 9, a signal 
envelope decision test 10, and two separate classification stages - one for processing 
signals with a constant envelope (11 and 13) and the other for processing signals with an 
irregular envelope (12 and 14). The purpose of the pre-classification stage 9 is to (a) 
determine if sufficient signal power is present, (b) estimate the signal's bandwidth and 
cenfr aid frequency, (c) perform the first level of carrier frequency error correction, and (d) 
filter the out-of-band background noise power. The first classification test is performed in 
decision step 10 where constant envelope signals are discriminated from those of irregular 
shapes. If the signal has a constant envelope, its modulation format is categorized within 
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the set {CW, FSK, FM}. Conversely, if significant envelope variations are detected, the 
signal's modulation format is categorized within the set {AM> DSB-SQ BPSK, QPSK, 
7T/4-QPSK, MPSK* OTHER}. Following this test, the residual carrier frequency errors 
superimposed on either signal are cancelled when required, in the error-correction steps 11 
and 12 respectively, and the processes far examining, recognizing and classifying the 
modulation schemes on the respective constant envelope or irregular envelope signals are 
performed in the classifiers 13 and 14, The phase and frequency content of the constant 
envelope signal is examined in classifier 13 in order to classify one of the {CW, FSK, 
FM} modulation formats. The phase and amplitude content of the irregular envelope 
signal is examined in classifier 14 in order to classify one of the {AM, DSB-SC, BPSK, 
QPSK, n/4-QPSK, MPSK, OTHER} modulation formats. 

Pre-classiflcation Stage 

Figure 2 shows in a flow chart a more detailed view of the Modulation 
Classification Section 3. The process begins at input IS where a block of digital complex 
samples is obtained from the output S of the Receiver Section of Figure 1. These samples 
are forwarded to a signal detector 16 where Fast Fourier processing determines if 
sufficient signal power is available to establish the presence of a signal. This process is 
performed with the aid of an external noise power module 17 that establishes a power 
threshold against which the power spectral density of the observed signal is compared. 
The noise power 17 is also used in by the B\V and centroid estimation process 20 to set 
another threshold for computing the centroid frequency and estimating the bandwidth of 
trie signal As illustrated in Figure 3, the noise power 17 corresponds to the level P w . A 
signal is declared present at decision step 18 if it bears some frequency components above 
a threshold (P^ ■+■ Z) dB, where Z is a given value relating to the probability that a noise 
frequency component has a power level above the (P w -r- Z) dB threshold. Thus, if the 
signal energy is not above the noise-relaxed threshold, the SNR is considered too small to 
allow further processing and the signal is classified as NOISE 19. However* if a signal is 
declared present, then the gross estimate of the carrier frequency inherent in the signal is 
determined in the B W and centroid estimation step 20. This operation is done by 
computing the centroid of the portion of the power spectral density that is above a certain 
threshold (Pw ■+■ Y) dB, as shown in Figure 3. This threshold is lower than the one used in 
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step 18 to detect the signal presence. Using Fast Fourier processing, the centroid frequency 
is given by: 

■ PSD[i] 

f =^ • < X > 



where PSD[# is the power spectral density value corresponding to frequency fi, and the 
summation is computed for the frequencies corresponding to power spectral density values 
PSD[f] above the threshold (P*, + Y)dB. Frequency translation of the signal is then 
performed in the gross error correction process 21 where the carrier frequency is corrected 
by an amount equal to the centroid frequency. The frequency bandwidth is also estimated 
in step 20. This operation is done by selecting the bandwidth corresponding to a pre- 
selected percentage of the sura of the power spectral density values PSD(i] above the 
given threshold (P w + Y)dB. This estimated bandwidth is then used in step 22 to filter the 
out-of-band noise power from the frequency-translated signal in step 20. The filter is 
selected from a bank of pre-stored filters. 

The process for detecting aod classifying the modulation format embedded in the 
signal it performed in steps 23, 30 and 31. Rapid and precise classification is possibLe 
provided that the residual carrier frequency error inherent on the output signals 28 and 29 
is kept to withm 0.001% of the sampling frequency F s used in sampling step 6 of Figure L 
This limitation is illustrated in Figure 4, where the success rate for four different 
modulation types is maintained at the zero-offset levels for frequency errors up to 10" s x 
F fi . It should be noted, however, that the presence of amplitude variation on the complex 
baseband signal is insensitive to any phase vector rotation caused by the inherent carrier 
frequency error. The present invention takes advantage of this feature when discrariinatmg 
between constant envelope signals and irregular envelope signals in decision test 23. This 
binary classification test is performed by using the maximum of the squared Fourier 
Transform on the normalized signal amplitude, given by 



|DFT( a )l- m 
r x = rnax — — ~ W 
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where/ is the frequency, DFT(.) is the Discrete Fourier Transform (computed from the 
Fast Fourier Transform), N s is the number of samples in the input block, and a is the 
amplitude vector centred on zero and normalized by its mean. Mathematically, vector a is 
expressed as 

where x is the observed filtered vector, x is the vector of magnitudes for the elements of 
vector x, and E[|xj] is the average of the magnitude elements over the vector jl The test 
parameter y m ^ is totally independent of the phase variation in x and is therefore 
unaffected by residual carrier .frequency errors. At the same time, it allows for the 
identification of information in the signal envelope, which, in most cases of interest, 
results in the reliable discrimination between constant envelope and irregular envelope 
signals. The performance obtained in Equation 2 is very good for SNRs as low as 5 dB, 
where the SNR is defined over the sampling bandwidth, 

Frequency Error Correction of Constant Amplitude Signals 

When test 23 determines that the signal has a constant envelope, the process moves 
to the error-correction step 26 where the residual carrier frequency error is estimated and 
corrected, The purpose of step 26 is to cancel out the carrier frequency error inherent in 
the constant amplitude signal 24 in order to determine if the signal is a continuous wave 
(CW or pure tone). Recognition and classification of a. CW signal cannot be performed on 
the constant envelope signal unless its inherent carrier frequency error has been corrected. 
The operation is performed in the frequency domain as shown in Figure 5, consisting of 
three signal proce$$ing operations: a coarse estimation of the earner frequency ( 32, 33, 34 
and 35); a refined estimation of the carrier frequency ( 37); and the carrier frequency error 
correction of the con$tant envelope signal ( 39 and 40). The process begins at step 32 
where the normalized squared signal is computed. The result is then processed by a Fast 
Fourier Transform (FFT) 33 with zero padding, squared in step 34, and forwarded to a 
coarse frequency estimator 35 where a search is done to locate the frequency line with 
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maximum power. Note that both the zero padding and the squaring of the output of the 
FFT improve the frequency resolution of the coarse frequency estimator 35, The output 
36, denoted by /max, is forwarded to a fine frequency estimator 37 where the refined 
estimation of the carrier frequency is performed. This signal processor (37) also receives 
the constant envelope signal 24 and computes the estimated residual carrier frequency 
error Af at 38, as follows: 

/ • PSPDL/^ ] - f m • PSDt>[/ ] (3) 
1 PSDDt/^l-PSDDE/J 



I? 5 where FSDD[/] is the derivative of the power spectral density at frequency/, and f x is a 
10 frequency given by 

fa fi-Ux±FJ2N 

p5 where F$ is the sampling frequency and AT the number of samples used in the FFT (that is, 
including zero padding). The plus sign is selected for// if PSDDI/naxl is positive, and the 
1 5 negative sign is selected if PSDD[£nax] is negative, PSDD[.] is obtained as 



PSDDLfl - 2Re[^ (/) J* (/)] 



^ where Re[_] indicates the real part, S(f) is the Discrete Fourier Transform (DFT) of the 
20 filtered observed signal at 24, S*(f) is its complex conjugate^ and S r (J) is its first derivative 
with respect to frequency, given by 



AT 



Jb=i 



2$ where x k is a sample, of tbe filtered observed signal vector at 24 and N is the nuiriber.of 



f 



samples in the input vector. 

The estimated frequency error Af at 38 is used to set the frequency of a digital 
VCO 39. The output of this VCO is then multiplied by the constant envelope signal 24 in 
a multiplier 40 to cancel out the signal's residual carrier frequency error. 
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Classification of Constant Envelope Signals 

The output 28 of Figure 5 (and Figure 2) is a frequency-translated, constant 
envelope signal. This signal is now forwarded to classifier 30 of Figure 2 where its 
embedded modulation is further examined to determine if it is CW, FM or FSK. This 
classification process is illustrated in Figure 6 and consists of two steps: (a) a test for 
significant phase variation to determine if the signal is CW or frequency modulated ( 41), 
and (b) a series of tests for deterrnining specific frequency characteristics to enable 
discrimination between FM and FSK formats (42, 43, 44, and 46). The process begins at 
decision step 41 where the frequency-translated, constant envelope signal 28 is examined 
to determine if significant phase variations are present within its envelope. Since CW 
signals have a near-zero variance in the unwrapped phase (direct phase) and frequency- 
modulated signals have relatively high variances, a test that detemrines the variability of 
the unwrapped phase is used. However, this unwrapped phase test requires care when the 
signal vector amplitude is small because phase variations increase as the SNR decreases. 
The difficulty of ibis problem is mitigated by discarding any data sample whose amplitude 
is below a threshold. Such threshold lowers the probability of residual phase variance in 
the CW signal, while noc affecting the phase variance of frequency-modulated signals. Ad 
amplitude threshold equal to the mean of the amplitudes of the sampled signal is chosen 
under the preferred embodiment of this invention. Discrimination between CW and 
frequency modulated signals js obtained by first, retaining the samples of the signal 28 at 
an amplitude above the amplitude threshold, and then comparing the value of the 
unwrapped phase of these samples to a phase threshold whose value is generally optimized 
through computer simulation experiments. If little or no phase variation exists, decision 
test 41 determines that the signal is CW and notifies the phase processor 42 not to proceed 
with further classification. Conversely, if significant phase variations occur in the signal, 
step 41 declares the presence of frequency modulation and notifies the phase processor 42 
to proceed with the next stage of classification process for discrimination between FM or 
FSK formats. 

To perform the classification between FM and FSK formats, the frequency- 
uncorrected, constant envelope signal 24 is used instead of the frequency-translated signal 
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28 because the process for computing the latter signal (described in Figure 5) may create 
large frequency errors for non-CW signals, and such large errors would create problems in 
determining the FM and FSK classifications. On the other hand, the residual carrier 
frequency error inherent in signal 24 is sufficiently small to allow for accurate 
discrimination between FM and FSK signals. In order to perform such discrimination 
between FM and FSK, the kurtosis coefficient of the instantaneous frequency content of 
the signal is used. This coefficient represents a measure of the flatness of the signal 
frequency distribution, and its value for analog FM is usually different from that of digital 
frequency modulated signals. The kurtosis coefficient also has a good behaviour at low 
SNRs. The coefficient is given by: 

where^ft) is the instantaneous frequency (about the mean) at time and B[fft)] is the time 
average o£f(t) computed over the length of the observed signal. 

The process begins at the phase processor 42 where the instantaneous frequency is 
obtained by computing the phase derivative of the constant envelope signal 24. This 
computation is further illustrated in Figure 7, consisting of signal processing steps 49, 50, 
51, 52, and 53. The phase of the constant envelope signal 24 is computed in step 49 to 
produce the instantaneous phase signal 50, The instantaneous phase signal is then 
processed through a low-pass filter (51 and 52) in order to limit the noise enhancement 
effects of the derivative function in step 53. A simple computation in BW estimation 51 
roughly estimates the effective bandwidth of the phase signal which, in turn, selects from a 
library one of a small number of low-pass filters that has a cut-off frequency slightly 
higher than the signal bandwidth. The output of the derivative function 53 is the estimate 
of the instantaneous frequency 54. This estimated value is required in the computation of 
the kurtosis coefficient in step 44 of Figure 6. It is also required to detect the possible 
presence of a low-modulated commercial FM signal in decision step 43. The FM 
modulation of such a signal contains the presence of silent periods, and, when 
demodulated, it produces a single discrete tone. Other forms of FM signals also contain 
such tones, such as the analog AMPS signals. Thus, the test performed on the constant 
envelope signal 24 in step 43 serves the purpose of detecting the presence of such a tone. 
If the tone is present, the signal is declared to be FM (commercial or AMPS). Conversely, 
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if such tone is not present, the signal is cither narrowband FM or FSK Following this 
process, the kurtosis coefficient is tested. FM signals usually produce a kurtosis coefficient 
that is larger than 2.5. Conversely, the coefficient produced by FSK signals is smaller 
because such signals are more compactly distributed. Thus, if the kurtosis coefficient is 
below the 2.5 threshold, a decision is made in favour of FSK. When the kurtosis 
coefficient is higher than the 2.5 threshold, an initial decision is made at 45 in favour of 
FM. However, a further classification test is required at this point because FSK signals " 
with highly filtered data (such as GMSK signals) may also have a high kurtosis 
coefficient. A frequency-modulated signal whose digital modulation contains a symbol 
rate off sym and a frequency deviation of a multiple of/^™/4, produces discrete frequency 
lines in the power spectral density (FFT) gf the squared signal. For example, MSK and 
GMSK signals have a frequency deviation of 2 x £ym/4 -f sym 12 and a modulation level 
M = 2. The FFT of the square of either signal produces two discrete frequency lines. 
Accordingly, the discrimination process is further refined in decision test 46 where the 
FFT of the squared constant envelope signal 24 is computed. (This computation is the 
same as the one done in steps 33 and 34 of Figure 5). Thus, if two or more discrete 
frequency lines are delected, the signal is classified as FSK at 48. Conversely, if decision 
test 46 produces less than two lines, the signal is classified as FM. 

This concludes the detailed description for the classification tests performed on 
constant envelope signals. 

Frequency Error Correction of Irregular Envelope Signals 

The following steps describe the classification processes for signals that have 
irregular envelopes. When decision step 23 of Figure 2 determines that the measured 
signal has an irregular envelope, the process moves to the error-correction processor 27 
where the residual carrier frequency errors are estimated andj corrected. This operation is 
performed in the frequency domain as illustrated irt Figure Sj and is similar to the 
operation described in Figure 5. Again, the process consists of three signal processing 
steps: a coarse estimation of the carrier frequency (55, 57, 59> and 6f); a refined 
estimation of the carrier frequency (62); and the carrier frequency error correction of the 
irregular envelope signal 25 (65 and 66). The process begins at step 55 where the square of 
the irregular envelope signal 25 is performed. The resulting signal is then processed by 
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the FFT 57 with zero padding, squared in step 59, and forwarded to a coarse frequency 
estimator 61 where a search is done to locate the frequency line with maximum power. 
Zero padding and the squaring of the EFT output 58 perform the same resolution functions 
as those described in Figure 5. Depending on the modulation format of the unknown 
5 signal, a peak in 'the frequency will appear in path 60. Choosing the frequency line with 
the highest peak, the process moves to a fine frequency estimator 62, where the carrier 
frequency estimation is further refined- This signal processor (62) also receives the 
normalized squared version 56 of the irregular envelope signal 25 and, using the one-step 
of the secant method as described in equation 3, computes the estimated carrier frequency 
L0 error. This value is halved in step 63 according to the squaring non-linearity in step 55. 
The result Af at 64 is used to set the frequency of a digital VCO 65 whose output {exp(- 

is then multiplied by the irregular envelope signal 25 in a multiplier 66 to cancel 
out the signal's residual carrier frequency error. As illustrated in Figure 9, the foregoing 
frequency error correction process is capable of obtaining a normalized RMS frequency 
15 estimation error lower than 10" 5 for the SNRs of interest (ie., larger than 5 dB). This carrier 
recovery method is therefore sufficiently precise to avoid performance degradation due to 
the limitations described in Figure 4. However, it should be noted that the frequency error 
estimation scheme performed in Figure 8 does not produce an adequate estimation of the 
earner frequency for MPSK signals. This limitation notwithstanding, the follow-on 
20 process leading to the classification of MPSK uses the irregular envelope signal with 
residual carrier frequency errors at 25 of Figure 2. 

Classification of Irregular Envelope Signals 

P The frequency-translated, irregular envelope Signal 29 of Figure 8 (and Figure 2) is now 

forwarded to classifier 31 of Figure 2 where it is further examined to determine if it is amongst the 
25 set (AM, DSB-SC, BPSK, QFSK, MPSK. OTHER} - This classification process is illustrated in 
Figure 10 and comprises the following signal processing steps: discrimination between one- 
dimensional and two-dimensional signal* (67): classification of one-dimenSional signals (68 and 
70); and classification of two-dimensional signals (73, 75, and 77). The process begins at decision 
test 67 where the frequency-translated, irregular envelope signal 29 is examined to determine if 
30 any modulation information is carried in the phase of the signal. If such information is not present, 
the signal modulation corresponds to the actual baseband, such as AM (transmitted carrier), DSB- 
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SC. or BPSK, and is recognized as one-dimensional. These signals are recognized from their 
absolute centred unwrapped phase sequence which has a small variance, Ir. practice, the process 
requires phase unwrapping because the initial carrier phase is random. However, phase 
unwrapping is very sensitive to noise, particularly on DSB-SC and BPSK signals due to theix 1 80 p 
5 phase transitions. Moreover, any phase unwrapping error will seriously increase the variance of the 
absolute centred phase. To avoid phase unwrapping of DSB-SC and BPSK signals, a new quantity 
is introduced in step 67 called the "absolute" phase, and is defined as 

A(0Ml/(0|+;|fl(0|) 

^ 10 where I(t) and Q(t) are the inphase and quadrature samples at time r. By taking the 

absolute values of the real and imaginary parts of a signal, a phase sequence between 0 
and til radians is produced. The phase sequence allows for a significant reduction in the 
size of the observation space required to make a decision. For DSB-SC signals the 
variance of the phase of the resulting signal is generally small. Small variance is also true 
15 for BPSK and AM signals. For two-dimensional signals bearing some phase infonnation 
(such as QPSK and MPSK signals), the variance of the absolute phase is high, tending 
towards the variance of a uniformly distributed random variable in the interval [0, rf2] 
radians.. To provide a better separation between one and two-dimensional baseband 
signals, a threshold on the amplitude of the signal is used in step 67, and the samples 
r u 20 below this threshold are discarded. The use of the amplitude threshold reduces the 

variance of the phase on one-dimensional signals, without significantly affecting the 
variance on two-<limensional signals. Simulations show that the best results occur when 
1 3 the sa id threshold is set equal to the mean of the amplitude of the phase sequence. If 

decision step 67 determines that the signal modulation is one-dimensional, the process 
25 moves to decision test 68 where the variance of the unwrapped phase (direct phase) is 

computed. Due to the presence of n radian jumps in the instantaneous phase of DSB-SC 
and BPSK signals, the said variance is much higher for these signals than for AM signals. 
Phase unwrapping is useful in this case because errors of this quantity are unlikely for AM 
signals, thus providing a good discriminating feature. In order to reduce the phase variance 
30 for AM signals, a threshold equal to the mean amplitude is set on the amplitude of the 

signal samples in step 68. An indication of low phase variance thus classifies the signal as 
AM at 69. Note also that ASK modulation is a digital form of amplitude modulation. ASK 
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signals arc therefore classified as AM signals at 69. If the phase variance is high, decision 
step 68 determines that the signal is not AM and moves to decision step 70 where a further 
test is made to determine if the signal is BPSK or DSB-SC. Unlike BPSK signals, DSB- 
SC signals bear an envelope for which the amplitude varies substantially over time. 
Accordingly, step 70 computes the variance of the envelope and uses it to discriminate 
between a DSB-SC signal at 72 and a BPSK signal at 71. 

Returning to step 67, if the test determines that the signal is two-dimensional the 
process moves to decision steps 73, 75 and 77 where PSK signals are separated from 
QAM and other unidentified modulation types. To perform the classification between 
QPSK, a/4-QPSK, MPSK and OTHER types of two-dimensional formats, the frequency- 
uncorrected, constant envelope signal 25 from decision step 23 is used instead of the 
frequency-translated signaL 29 computed in Figure 8, because the residual carrier 
frequency error inherent in signal 25 is sufficiently small to allow for accurate 
discrimination between the QPSK, tt/4-QPSK, MPSK and OTHER signals. To initiate this 
classification in decision step 73 7 a simple test on the variance of the signal amplitude is 
performed. If the variance of the amplitude is below a given threshold, the signal is 
assumed to be a PSK signal. Even though band-limited PSK signals do exhibit amplitude 
variations, such variations are much less noticeable than those for QAM, SSB or VSB 
signals. This test is therefore sufficient to discriminate PSK signals from most of the other 
two-dimensional signals as OTHER at 74, If the test determines thai the signal is PSK, the 
process moves to decision step 75 where the signal is further examined to determine if its 
embedded modulation is QPSK, n/4-QPSK or MPSK. This test is performed by 
computing the fourth power of the signal 25. QPSK signals produce a single peak in the 
power spectral density of the resulting signal, while tt/4-QPSK signals produce two peaks 
separated in frequency by twice the baud rate. Therefore, if the test in step 75 results in 
one peak,, the signal is classified as QPSK at 76. Otherwise the process moves to decision 
test 77 where the combination of two peaks classifies the modulation as tc/4-QPSK at 78 
and the absence of any peak classifies the modulation as MPSK at 79. 

This concludes the detailed description for the classification tests of irregular 
envelope signals. The description covering the overall decision tree process as illustrated 
in Figure U , is given below. 
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O verall Description of the Preferred Embodiment 

Figure 1 1 illustrates a complete overview of the modulation classification section 3 
of Fieure 1 and collects together the signal processing functions illustrated in Figures 2 ? 6 
and 10. The process begins at input 15 where a block of digital complex samples is 
obtained from the output 8 of the Receiver Section of Figure 1. These samples arc 
forwarded to a signal detector 16 where Fast Fourier processing determines if sufficient 
signal power is available to establish trie presence of a signal. This process is performed 
with the aid of an external noise power module 17 that establishes a power threshold 
against which the power spectral density of the observed signal is compared. The noise 
power module 17 is also used in step 20 to set another threshold for computing the 
centroid frequency and estimating the signal bandwidth of the signal. If the signal energy 
is not above a noise-related threshold, the SNR is considered too small to allow further 
processing and the signal is classified as NOISE in step 19, However, if a signal is 
declared present, then the gross estimate of the carrier frequency inherent in the signal is 
determined in step 20. This operation is done by computing the centroid of the portion of 
the power spectral density that is above a certain threshold. Frequency translation of the 
signal is then performed in correction step 21 where the carrier frequency is corrected by 
an amount equal to the centroid frequency. The frequency bandwidth is also estimated in 
estimation step 20. This estimated bandwidth is then used in filter step 22 to filter the out- 
of-band noise power from the frequency-translated signal in estimation step 20. The 
process for detecting and classifying the modulation format embedded in the signal starts 
in decision test 23 where discrimination between constant envelope signals and irregular 
envelope signals is performed- When decision test 23 determines that the signal has a 
constant envelope, the process moves to the error-correction step 26 where the residual 
carrier frequency error is estimated and corrected. The classification of constant-envelope 
signals starts in decision test 41 and consists of the following steps: (a) a test for 
significant phase variation to determine if the signal is CW or frequency modulated (41), 
and (b) a series of tests fox determining specific frequency characteristics to enable 
discrimination between FM and FSK formats (42, 43 7 44 and 46). If little or no phase 
variation exists in the frequency-translated signal, decision test 41 determines that the 
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signal is CW and notifies the phase processor 42 not to proceed with further classification. 
Conversely, if significant phase variations occur in the signal, decision test 41 declares the 
presence of frequency modulation and notifies the phase processor 42 to proceed with the 
next stage of classification process for discrimination between FM or FSK formats. To 
5 perform the classification between FM and FSK formats, the frequency-uncorrected, 
constant-envelope signal 24 from decision step 23 is used. The process begins at 
computing step 42 where the instantaneous frequency is obtained by computing the phase 
derivative of the constant envelope signal 24. This estimated value is required in the 
computation of the kurtosis coefficient in decision test 44. It is also required to detect the 
10 possible presence of a low-modulated commercial FM signal in decision test 43. The FM 
modulation of such a signal contains the presence of silent periods, and, when 
demodulated, it produces a single discrete tone. Other forms of FM signals also contain 
such tones, such as the analog AMPS signals. Thus, the test performed on the constant 
Hj envelope signal 24 in decision test 43 serves the purpose of detecting the presence of such 

0 15 a tone. If the tone is present,, the signal is declared to he FM (commercial or AMPS). 

Conversely, if such tone is not present, the signal is either narrowband FM or FSK. 
Following this process, the kurtosis coefficient is tested in decision test 44. FM signals 
usually produce a kurtosis coefficient that is larger than 2.5. Conversely^ the coefficient 
;<«i produced by FSK signals is smaller because such signals are more compactly distributed- 

!' il 20 Thus, if the kurtosis coefficient is below the 2,5 threshold, a decision is made in favour of 

FSfL When the kurtosis coefficient is higher than the 2.5 threshold, an initial decision is 
H made in favour of FM. However, a further classification test is required at this point 

13 because FSK signals with highly filtered data (such as GMSK signals) may also have a 

high kurtosis coefficient The discrimination process is further refined fn decision step 46 
25 where the FFT of the squared constant envelope signal 24 is computed- If two or more 
discrete frequency lines are detected, the signal is classified as FSK. The estimate of the 
number of symbols M is obtained in counting step 47 by counting the number of FFT 
peaks. Conversely, if decision test 46 produces less than two lines, the signal is classified 
as FM. 

30 When decision step 23 determines that the measured signal has an irregular 

envelope, the process moves to the error-correction step 27 where the residual carrier 
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frequency errors are estimated and corrected. The frequency-translated, irregular envelope 
signal output of error correction step 27 is now forwarded to decision test 67 where it is 
further examined to determine if it is amongst the set {AM, DSB-SQ BPSK, QPSK, 
MPSK, OTHER). This classification process comprises of the following signal processing 
5 steps: discrimination between one-dimensional and two-dimensional signals (67); 
classification of one-dimensional signals (68 and 70); and classification of two- 
dimensional signals (73, 75, and 77). The process begins at decision test 67 where the 
frequency-translated, irregular envelope signal is examined to determine if any modulation 
information is carried in the phase of the signal. If such information is not present, the 

10 signal modulation corresponds to the actual baseband, such as AM (transmitted carrier), 
DSB-SC, or BPSK, and is recognized as one-dimensional. If decision test 67 determines 
that the signal modulation is one-dimensional ? the process moves to decision test 68 where 
the variance of the unwrapped phase (direct phase) is computed. An indication of low 
phase variance thus classifies the signal as AM. If the phase variance is high, decision test 

15 68 determines that the signal is not AM and moves to decision step 70 where a further test 
is made to determine if the signal is BPSK or DSB-SC. Unlike BPSK signals, DSB-SC 
signals bear an envelope for which the amplitude varies substantially over time. 
Accordingly, decision test 70 computes the variance of the envelope and uses it to 
discriminate between a DSB-SC signal from a BPSK signal. 

20 Returning to decision step 67, if the test determines that the signal is two- 

dimensional, the process moves to decision steps 73, 75 and 77 where PSK signals are 
separated from QAM and other unidentified modulation types. To perform the 
classification between QPSK, ti/4-QPSK, MPSK and OTHER types of two-<Jimensioiial 
formats, the frequency-uncorrected, constant envelope signal 25 from decision test 23 is 

.25 used. To initiate this classification in decision step 73, a simple test on the variance of the , 
signal amplitude is performed. If the variance of the amplitude is below a given threshold, 
the signal is assumed to be a PSK signal. Otherwise, it is classified as OTHER. If the test 
determines that the signal is PSK, the process moves to decision step 75 where the signal 
is further examined to determine if its embedded modulation is QPSK, tl/4-QPSK or 

30 MPSK This test is performed by computing the fourth power of the signal 25. QPSK 

signals produce a single peak in the pow er spectral density of the resulting signal, while 
Tt/4-QPSK signals produce two peaks separated in frequency by twice the baud rate. 
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Therefore, if the test in step 75 results in one peak, the signal is classified as QPSK, 
otherwise the process moves to decision step 77 where the combination of two peaks 
classifies the modulation as k/4-QPSK and the absence of any peak classifies the 

modulation as MPSK- 

Different aspects of the embodiment in Figure 1 are possible. For example, the 
signal coming to the Receiver Section 2 can be measured from a tapped wireline, or from 
another receiver or pre-amplifjcation section- The signal may also have been amplified and 
filtered by an alternate receiver, and be introduced before the frequency down conversion 

5. Similarly, the down conversion to complex baseband could be performed on an 
alternate set of hardware, and the signal could be introduced before the sampling function 

6. Furthermore, .the sampled digital complex baseband signal 8 could be obtained from a 
receiver section with a different configuration, or it could be retrieved from a digital 
memory location where it had been previously stored. Note that, for practical reasons, the 
sampling function 6 and the AfD converter 7 are usually implemented in a single unit, and 
that it is not very likely that the signal could be introduced before the digital conversion 
function 7. Such a procedure is nevertheless conceptually also possible. 

PERFORMANCE EVALUATION 

In order to demonstrate the methods described in the present invention, 500 
simulated signals of each of the modulation types have been generated and processed. 
These signals covered a wide variety of parameters as described in the next section, A 
sampling frequency of 48 kHz was used covering a bandwidth slightly larger than the 
occupied bandwidth of most narrowband communications signals. Sequences of 853 msec 
(representing 4096 samples) were used as inputs to the modulation classifier (block 3 of 
Figure 1). Complex baseband signals were used, whereas the carrier frequency and phase 
errors were simulated The frequency error was random and uniformly distributed over a 
range [-4.8 kHz to 4,8 kHz], while the initial carrier phase was uniformly distributed over 
[-7L to tt]. For digitally modulated signals, a random delay uniformly distributed over [0 to 
TsJ was used to simulate symbol timing uncertainties (where T s is the sampling interval 
equal to 1/48000 seconds). 
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Simulation of Specific Signals 

Analog Modulations 



For analog modulation schemes, two types of source signals were simulated. The first was 
a real voice signal, band- limited to [0 to 4 kHz], The second was a simulated voice signal 
that used a first-order autoregressive process of the form: 

// y[k]=0.95xy[k^l]+n[k] 

where n[k] is a white Gaussian noise process. Furthermore, this pseudo-voice signal was 
band4irnited from 300 to 4000 Hz. 

For AM signals, a constant value was added to the source signal. The modulation 
index was calculated by using the maximum amplitude value over the whole source signal. 
The index was then uniformly distributed in the interval [50% to 100%]. The total length 
of the real source signal was about 120 seconds, while that for the pseudo- voice was about 
40 seconds. From these two source signals, sequences of S5.3 msec in duration were 
randomly extracted. Thus a the observed modulation index for a sequence was equal to or 
less than the chosen modulation index. 

For frequency-modulated signals, a cumulative sum was used to approximate the integral 
of the signal source. Generic FM signals were simulated using real or pseudo -voice signals 
with a modulation index uniformly distributed in the interval [1 to 4], The bandwidth 
occupied by these signals ranged from 16 kHz to 40 kHz, using the approximation: 




where /? is the modulation index andj^x is the maximum source frequency (4 kHz in this 
case). The AMPS FM signals were approximated using a modulation index of 3. 

Digital Modulations 

Continuous-phase FSK signals were simulated by using filtered M-ary symbols to 
provide frequency modulation to a carrier frequency. Pager signal parameters were based 
on observations of real signals, with 2FSK modulation at a bit-rate of 2400 bps, a 
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frequency deviation of 4.8 kHz, and almost no filtering. 4FSK signals were also simulated, 
using the same 4.8 kHz frequency deviation and a symbol rate of 1200 baud. Tbe 19.2 
kbps, 2FSK signals from theRacal Jaguar radio were simulated, using a frequency 
deviation of 6.5 kHz and a 5* order Butterworth pre-modulalion filter with a cutoff 
frequency of 9.6 kHz. Also included in the simulation were GMSK signals that were 
similar to GSM signals having a BT product of 0-3. 

For PSK and QAM signals, the symbols were filtered with either a raised cosine 
function or a square root raised cosine function. The selection was randomly performed 
with equal probabilities. The rolloff factors of 20%, 25%, 30%, 35%, 40%, 45%, and 50% 
were uniformly and randomly selected. For all these signals, the symbol rates were chosen 
randomly between 16 and 20 kbaud. Also simulated were W4-QPSK signals that were 
similar to IS-54 signals, with a symbol rate of 24 kbaud and a square-root raised cosine 
pulse-shaping filter with a rolloff factor of 35%. 

Additive Noise 

The simulated signals were passed through an additive white Gaussian noise 
channel before being classified. For the simulation, no filtering was done at the receiver, 
therefore the signal observed by the modulation classifier was corrupted by the white 
noise. The noise power was calculated from the knowledge of the average power of the 
modulated signal and the SKR over the sampling bandwidth. This SNR was defined as: 



I SNR ^mp = S/(N 0 . FJ 

where S is the signal power, N 0 is the white noise, power spectral density, and F 5 is the 
sampling frequency equal to 48 kHz. 

i 

For amplitude-modulated signals [AM, DSB-SC, and SSB], the amplitude power 
was calculated by using all source signals (real and pseudo-voice). This condition implies 
that, for a given sequence, the observed SNR might be different from the overall SNR, 
which is especially true for DSB-SC and SSB signals where some segments of the signals, 
because of silence segments, may have no power at all. 
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Binary Decision Thresholds 

With respect to the decision tree analyses performed under the present invention, 
each decision compares a signal feature with a mreshold. For the simulations undertaken, 
the thresholds were set from the results of direct observations of the feature distributions 
in a training set of simulated signals having anSNR of 5 dB. These selected thresholds are 
summarized in Table 1 for the different features. 



f ■ — "~ ~ 

| Feature 


Threshold 


[power level above noise level, for detection of signal presence (Step 16) 


10 dB 


1 Power level above noise level, for bandwidth and centroid estimation 
| (Step 20) 


5 dB 


Maximum of the normalized squared FFT of the centered normalized 
envelope (dma*) (Step 23) 


1.44 


Variance of the direct phase (tad) (CW vs. FM/FSK in Step 41) 


0.25 


| Variance of the direct phase (rad) (AM vs. DSB-SC/BPSK in Step 68) 


4.0 


Variance of the " absolute" phase (rad) (one-dimensional vs cvvo- 
1 dimensional signals in Step 67) 


0.144 


1 Kurtosis of the instantaneous frequency (FM vs FSK in Step 44) 


2.5 


(Variance of the normalized amplitude in one-dimensional signals (DSB- 
SC vs BPSK in Step 70) 


0.25 


1 Variance of the normalized amplitude in MP SK signals (OTHER vs 
1 QPSK/ti/4-QPSK/jWPSK in Step 73) 


0.15 



Table 1: Decision thresholds used in the simulations. 



D Classification Results 

An estimate of the performance of the modulation classification method of the 
present invention was obtained by applying the simulated signals described earlier. For 
each modulation type, the 500 generated sequences were classified by using the preferred 
decision tree mediods illustrated in Figure 11. Eleven outputs from the modulation 
5 classifier were possible: NOISE, CW, AM, DSB-SC, FM, FSK, BPSK, QPSK, n/4-QPSK, 
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MPSK, and OTHER Table 2 shows the classification results for each of the simulated 
modulation types at a SNR of 5 dB. The results of this table were obtained in the presence 
of uniformly distributed frequency errors between - 4.8 kHz and 4.8 kHz (10% of the 
sampling frequency). The initial carrier phase was also uniformly distributed over the 
range of ^ to ti. For digitally-modulated signals, a random delay uniformly distributed 
over one sample period, was used to simulate symbol timing uncertainties. 
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6.0 
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1-6 
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52.6 


SSB 
(SV) 












! 










100,0 



Table 2 - Classification results (in %) for a SNR of 5 dB, 



In Table 2, it is important to note the difference between the classification for 
analog modulation signals using real voice (denoted by •» V) and those using simulated 
voice (denoted by "SV"). For these modulations, the real voice signals included pause and 
silent durations, whereas the simulated voice signals were generated continuously, 
without any silent duration. For analog modulation, the pauses in a real voice source 
produced unmodulated sequences. For DSB-SC and SSB modulated signals, such pauses 
produced no signal at all. In the case where the duration of the 85.3 msec sequence was 
mostly a pause, the signal was classified as noise or as OTHER. In the case where the 
transition between a pause and voice was not clear within the 85.3 msec observation, 
erroneous modulation types appeared (as the row for SSB(V) in table 2 indicates). For AM 
and FM signals, the absence of a source signal produced a CW classification. If the 
duration of die silence occupied most of the S5.3 msec observation time, such quiet 
sequences were classified as CW signals. Again, depending on the transition timn between 
a pause and a voice, strange results appeared (as the row for FM(V) of Table 2 indicates). 
These results are not considered classification errors as such. Rather, they reflect the fact 
that the classification is more difficult to perform for analog-modulated signals when 
using a very short observation time. Such perceived errors would be elimhiated by either 
increasing the length of the observation time before the classification i 5 undertaken, or 
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adding a post-processing step that accumulates the results of several observations and 
performs a decision according to the dominant modulation type. 

Modular Construction 

The modular nature of the classification process illustrated in Figure 1 1 allows for 
different modulation schemes to be classified. Furthermore, although the disclosed 
apprgach under the present invention assumes an AWGN channel, there are possibilities of 
extending the channel to more complex models by employing additional processing, such 
as the m etho d of blind equ alization . 

Changes and modifications in the specifically described embodiments can be 
carried out without departing from the scope of the invention which is intended to be 
limited only by the scope of the appended claims. 
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